Improved Superconducting Qubit Readout by Qubit-Induced Nonlinearities 
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In dispersive readout schemes, qubit-induced nonlinearity typically limits the measurement fidelity by 
reducing the signal-to-noise ratio (SNR) when the measurement power is increased. Contrary to seeing 
the nonlinearity as a problem, here we propose to use it to our advantage in a regime where it can increase 
the SNR. We show analytically that such a regime exists if the qubit has a many-level structure. We also 
show how this physics can account for the high-fidelity avalanchelike measurement recently reported by 
ReedefflZ. l arXiv:1004.4323i /ll. 
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Quantum measurements are crucial to all quantum infor- 
mation protocols. In solid-state systems, readout can be 
performed by connecting the qubits to noisy measurement 
electronics, such as single-electron transistors [jj]. Disper- 
sive readouts based on coupling qubits to high-Q resonators 
are however much less disruptive as all of the energy is dis- 
sipated away from the qubit f3\. This type of measurement 
leads to quantum nondemolition (QND) readout. Unfor- 
tunately, so far the typical signal-to-noise ratio (SNR) has 
been relatively low, requiring sometimes up to 10^ repeti- 
tions of the experiment to average out the noise added by 
the amplifiers yl |j]. Increasing the qubit-resonator cou- 
pling is one approach to improve the SNR as it increases 
the amount of information about the qubit state carried by 
the measurement photons. This is however at the cost of an 
increased Purcell decay rate [SJ], which decreases the use- 
ful integration time and in turn the SNR. Another approach 
has been to exploit bifurcation in a nonlinear resonator ml. 
This has already led to the experimental demonstration yji 
of single-shot readout of a transmon-type superconducting 
qubit 181]. Very recently. Reed et al. have also shown that 
single-shot measurement of a transmon qubit can also be 
realized in a linear resonator but working at very large mea- 
surement power ||9|]. 

Motivated by these results, we study how the nonlinear- 
ity induced by the qubit in a linear resonator can lead to 
improvement of the measurement. We first explore the low 
driving power regime before considering the high-power 
regime studied in Ref. ^. Using a simple model, we find 
qualitative agreement with these experimental results. 

For simplicity, we first focus on a two-level system 
(TLS), of states labeled { 1 0) , 1 1) }, dispersively coupled to 
a microwave resonator. This system is well described by 
the Jaynes-Cummings Hamiltonian, expressed here in the 
dispersive basis {h = 1) |1( 
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where Hij = \i) {j\ for the TLS, and N = a^a. In this 
expression, w,. is the resonator frequency and cj^ (CJa) the 
bare (Lamb-shifted) TLS transition frequency. The qubit- 
resonator coupling is characterized by x = 5^(1 — '^^)/A 
the dispersive coupling strength and ( = —gX"^ is a Kerr- 
type nonlinearity, with g the bare qubit-resonator coupling 
strength, A = oji^ — w,. the qubit-resonator detuning and 
A = —g/A. The second line of Eq. ([B is valid to fourth 
order in A and at photon numbers n <C ?^crit> where rid it = 
1/4A^ is the critical photon number pfl. The limit C ~^ 
of the second line of Eq. ([T]l is the standard linear dispersive 
Hamiltonian ifllll . 

In this regime, because of the qubit-dependent pull of the 
resonator frequency x^z, a nearly resonant drive on the res- 
onator will displace the resonator's initial vacuum state to 
a qubit-state dependent coherent state 1 0^=0,1)- Homodyne 
measurement of the transmitted or reflected signal can then 
resolve these pointer states, and hence the qubit states. In 
Ref. 11211 , it was shown that the SNR of such a homodyne 
measurement for an integration time Ti = I/71 is given by 
SNR = rjK\ai — 0:0^/71, where k is the resonator decay 
rate and r] is the measurement efficiency. The SNR reaches 
its maximal value (SNR),„ax = 4r/nx/7i, for the optimal 
choice K = 2x, where n is the average number of photons. 
In this limit, the dispersive model Eq. ([T]) thus predicts that 
the SNR should increase linearly with the number of mea- 
surement photons n. 

Unfortunately, this linear increase of the SNR is not ob- 
served experimentally and this can, at least partially, be ex- 
plained by the fourth order corrections in Eq. ([T]l. To this 
order of approximation, an increase of the photon number 
leads to a reduction of the cavity pull ffA[l — A'^(a^a-|-1)](7^ 
and hence to a loss of distinguishability of the pointer 
states llioll . Thus, the prospects for improving the SNR in 
homodyne measurement of a TLS dispersively coupled to 
a resonator look rather unpromising. 

Fortunately, most superconducting qubits are well ap- 
proximated by many-level systems (MLS), often with only 
weak anharmonicity ll8L [l3l4l5n . rather than by TLS. As is 
apparent below, it is possible in this situation for x and ^ to 



have the same sign, which yields an increase of the cavity 
pull with n and hence an improved SNR. 

As a good description of a generic superconducting 
qubit, we thus consider an M-level system, of states la- 
beled {|0),...,|M-1)} and with the first two states 
acting as logical states. The Hamiltonian of the MLS- 
resonator system takes the generalized Jaynes-Cummings 
form H, = Ho + ^,^0" gi{a'^Ili^,+i + aU.+i^i), where 

Hq = ojrCi^a + Yli=Q ^i^ii i^ the free Hamiltonian, Ui 
the frequency of level |i), and gi the coupling strength be- 
tween the resonator mode a and the i -H- i + 1 MLS transi- 
tion. 

In the dispersive regime, where 4(a^a)(gi/A,)^ <C 1 
with A,; = {uji^i — oJi) — 00,., it is convenient to approxi- 
mately diagonalize Hg. Following Ref. lUoll where this was 
done for a TLS, we obtain to fourth order in A^ = —gij^i 



M-l 
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iff ^H,+ Y^ 5,n,,,ata + ^ K.U^^a^a) 
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where Hq includes the Lamb shifts. In this expression, we 
have defined the ac-Stark (Si) and Kerr (Ki) coefficients 

S, = [x,-i(l - A^) - x.(l - XU) - 2x.-iXU] 

+ i(9Xi-2A,^_i - 3xt-iXf-2 - X»A-+i + 3x^+iA- ) 
-^fAf-3,EA!!^„ (3a) 

K^ = i(3xj-2A,^_i - x»-iA?_2 + XiXf+i - 3Xi+iA- ) 



+ te-X.-i)(A? + ALi)+5rA, 
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AiAi4-i(Ai+i — Ai), Aj 
gl 7(A,+i + A,), and x^ = X^ = Oiori(^ [0, M - 2]. 
For M = 2,So = -Xo, Si = Xo(l - 2Ag) and Ko = 
—Ki = XoAq, reproducing Eq. (dJ. The crucial feature 
of these expressions is that, contrary to a TLS, Kq 7^ —Ki 
for a MLS. Using this state-dependent nonlinearity, we now 
show how to improve the readout in two ways. 

Weak driving limit. From the dispersive model Eq. (O, 
it is possible to evaluate the difference in cavity pull 5 for 
levels |0) and 1 1) which we are interested in distinguishing 
in a measurement. It takes the form (5 = x' + C'^' with 
n = (^a'^a'j and where we have defined x' = Si — Sq and 
(' = Ki — Kq. We note that, although we are focusing 
on the first two levels, the presence of higher MLS levels 
is important. This is apparent in the expressions for Si and 
Ki which involve states up to \i + 2). 

Figured] shows x' ^^^ C obtained from Eq. dS]) as a func- 
tion of ujr for a MLS with (a) M = 2 and (b) M = 6 
levels. Figure 1(c) has been obtained from exact diagonal- 
ization of H, for M = 6. For M = 2, sgn(x')=sgn(C') 
only in regions where the dispersive approximation breaks 
down while for M > 2 this is possible in the dispersive 
regime (see the caption of Fig. 1). The numerical results 
illustrate that the analytical expressions are good approx- 
imations. We note that, while these results apply to any 
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FIG. 1: (color online) Analytical (a,b) and numerical (c) ac-Stark 
x' (black line) and Kerr shifts C,' [lighter (red) line] for a trans- 
mon qubit taking into account M = 2 (a) and M — & levels 
(b),(c). The parameters are chosen such that (a;io, LiJ2i,ga}/2-K = 
(6000, 5750, 100) MHz. The lines change from dotted to full 
when x' or C,' pass from negative to positive values. Vertical 
blue dotted lines indicate transitions wio, 1^21, and LJ32. Lighter 
(green) shading indicates regions of interest, while darker (red) 
indicates regions where the dispersive model breaks down even 
at low-photon number. Values of x' and C^' obtained numerically 
are not plotted close to divergences. 
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FIG. 2: (color online) Cavity pull 5 (a) and SNR (b)-(d), as a 
function of the average number of photons. The parameters are 
the same as in Fig [T] with T\ — \ fj,s. Dashed blue lines cor- 
respond to the second order approximation, for which the cav- 
ity pull 2g^/A is constant. Full black lines (dotted red lines) 
are obtained for ujr/2iT = 4515(7660) MHz, corresponding to 
"^crit ~ 55(70) and sgn(x') = {—)sgn{('). These frequencies are 
chosen such that |x'I/2t = 2 MHz. 



MLS, for concreteness we have chosen here parameters 
corresponding to a transmon qubit [SI]. 

Figure 13 a) shows the cavity pull 6 for a cavity fre- 
quency chosen in the region where sgn(x')=sgn((^') (full 
black line) and outside of this optimal area (dotted red line). 
These results are compared to the result of the second order 
approximation (dashed blue line). As expected, the cavity 
pull increases with n under the appropriate choice of pa- 
rameters. Figures [2tb)-|2d) show the corresponding SNR, 



using the same color scheme, for k,/2x' = 0.5 Eb)], 1 
lEtc)] and 1.5 ||2d)], with k/2x' = 1 being the optimal 
choice 11211 . Because of the increase of the cavity pull with 
n, the full black line is always above the dotted red one. For 
k/2x' = 1, an improvement of nearly 100% is expected 
at large photon numbers. Biasing the qubit above the res- 
onator's fundamental frequency, as is suggested here, can 
lead to an increase of the Purcell decay. This can however 
be strongly reduced by a small change of design flol . Fi- 
nally, one could also tune the system to a point where x' 
and Q' have the same sign at the moment of measurement 
using a tunable resonator lIlTlllSll . 

Strong driving limit. The results obtained so far relied on 
Eq. (O, which is valid only below ricrit. However, the non- 
trivial state dependence of the nonlinearity Ki should ex- 
tend well beyond the dispersive regime. To explore this, we 
take advantage of the block diagonal structure of the Hamil- 
tonian Hs- There each block of Hg corresponds to a fixed 
number n of qubit-resonator excitations and is spanned by 
£^ = {|n, 0) , ..., \n - M + 1, Af - 1)}. With the rele- 
vant M being at most ~ 10 in practice, we can diagonalize 
each block numerically (or analytically for M < 4) for ar- 
bitrary n. In this way, we obtain the dressed energies £'„_j 
and states In, i), where E„i is the energy of the eigenstate 
closest to a Fock state with n photons and MLS state \i). 
From these expressions, we find the effective resonator fre- 
quency uJriiiT-) = En+i,i — En.i- This frcqucncy depends 
in a nonlinear way on the MLS state-dependent average 
photon number iii. In steady state and in the absence of 
qubit transitions, rii is given by the measurement drive am- 
plitude e and frequency detuning relative to the effective 
resonator frequency 

n,; (e , W™ ) = -^—^ 1, , r /»in ' (4) 



\OJr 



with Lo^n the measurement frequency. We solve Eq. ^ iter- 
atively to find n^ and oj^i as a function of e and u}„i. 

In Fig. [21 both ujri l|3la)-13c)] and n, ||ld)-[3l;f)] are plot- 
ted as a function of measurement power for M = 2, 3, 6 re- 
spectively. This is done for w^ = Wr as in Ref. \9\. For all 
values of M, the effective resonator frequency approaches 
its bare value w,. at large power This is expected because, 
at this point, {N) » (V^) in Eq. ([T]l and the cavity re- 
sponds classically IJI^]. Since Kq = — A'l for M = 2, 
this classical crossover occurs at the same input power for 
both qubit states. We note that this crossover happens in 
an avalanche manner, with each additional photon bringing 
ujri{ni) closer to uj^ and facilitating the addition of more 
photons. As can be seen in Fig. [3](b,c), for M = 3 this 
avalanche occurs at a state-dependent power, the behavior 
changing only quantitatively for M > 3. With oj„i = u^, 
we thus expect an abrupt change in the average photon 
number in the resonator at a power that is MLS state de- 
pendent. This is illustrated in Fig. [S^f) where for M = 6 
there is a range of ~ 5 dB in measurement power where 
Uq and Til differ significantly and by as much as ~ 10^ at 
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FIG. 3: Effective resonator frequency uj^i (a)-(c) and mean pho- 
ton number n^ (d)-(f) for i = (full red lines), i = 1 (dotted blue 
lines), and i = 2 [dashed gray lines, (c),(f)] as a function of the 
measurement power. Panels (a,d), (b,e) and (c,f) are for M — 2, 
3, and 6 respectively. In (a)-(c), the dashed green horizontal line 
is ujr- The parameters are the same as in Fig[T] For clarity of 
presentation the panels have different horizontal scales. 
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FIG. 4: Mean photon number rii for i = (a) and i = 1 (b) 
as a function of the measurement frequency and power The full 
white lines are ujriirn) for Umjl-K — LJr/2-K = 7 GHz. The 
solid red and dotted blue vertical lines indicate the measurement 
frequency used in Fig[3l The parameters are the same as in Fig[T] 
The horizontal dotted lines delimitate the regime of measurement 
power where ni — jiq is maximum. 



the optimal driving power. This large separation of the S 
curves, much larger than typical amplifier noise, leads to 
single-shot readout of the qubit 131. As shown by the gray 
dashed line in Fig [51 f), pumping the 10 2 transition be- 
fore readout 13, 1^1 could also help the measurement work 
at lower power. Finally, Fig.|4]shows the full power versus 
measurement frequency dependence of iii (e, cj,„) obtained 
from Eq. @. In this plot, the full white lines correspond 
to ujri shown in Fig. Oc). Although Fig. 0] shows qualita- 
tive agreement with the results of Ref. ||9il, a quantitative 
comparison would require proper modeling of the 4-qubit 
device used in Ref. H. 

While the dispersive measurement at low-photon num- 
ber of Fig. |2]is expected to be QND, this is not the case 
for the single- shot high-power measurement. To eval- 
uate the QND character of this avalanche readout, we 



estimate how the presence of the measurement photons 
changes the relaxation and excitation rates of the qubit, 
as well as causes leakage outside of the logical subspace 
{|0),|1)}. Figure |5ja) shows the Purcell decay rate 

7^/k ~ I (n, Ol a In, l) I and corresponding leakage rate 
Yli=to 1 K'^ ~ ^'^1 '^ l'^' l)| ^^ a function of mea- 



li/^ 



surement power. At low power, we find the expected result 
7k /^ = -^0 iSl- For large photon number, n + 1 w n 
and losing a photon through resonator decay does not sig- 
nificantly change the qubit states. The Purcell decay rate 
thus goes down with measurement power and does not af- 
fect the QND character. To evaluate how the qubit dress- 
ing changes pure relaxation. Fig |3b) shows the rates for 

dressed decay 71^/71 ~ | (n, 0| S_ |n, l) | and leakage 



7id/7i 



n 



i, i\ $]_ n, 1 



where S_ 



M-2 



Eivi — 
i=0 



ffo 



Hi i+i. The participation of the higher transmon 



states reduces the decay rate 71^ from |1) to |0) as mea- 
surement power is increased. However, decay of the bare 
higher states increases the leakage rate, and the total error 
rate 7id + 7^^ is larger than 71. 

Finally, dressed-dephasing 7^ due to noise responsible 
for dephasing of the bare qubit states can also cause tran- 
sitions between the dressed-states lUol llOti . For concrete- 
ness, we consider dephasing due to charge noise on a trans- 
mon, but the model can be adapted to any source of dephas- 
ing. Following Ref. llOfl, this contribution can be evaluated 

as 7d/7^ « K^^+Ml S, |M>r S(Aio)/5(l Hz), 
where S^ = ^^^q Hjiej/ei, with e^ the charge dis- 
persion of level i ^ and S{/S.ij) the spectrum of 
charge noise evaluated at the dressed qubit-resonator de- 
tuning. In the same way, the leakage rate is 7^/71^ ~ 

X:,^o_, I {n + l-i,i\ S, |M> r S{Ku)/S{l Hz). 

We note that, even assuming 1// charge noise which 
would be 10*^ times smaller at 1 GHz than at 1 Hz, dressed 
dephasing can be important for these large photon numbers. 
Indeed, for the transmon, the charge dispersion e.i — and 
therefore the susceptibility to charge noise — increases ex- 
ponentially with i, reaching eg/ei ~ 10^ for 6 levels. Al- 
though a quantitative analysis requires a better understand- 
ing of the noise spectrum at microwave frequencies, our 
numerical analysis with 1// noise suggests ^d/l^ ranging 
from 1 to 10"^ and 7^/7,^ from 10 to 10'^ depending on the 
parameters and the number of levels. 

The loss of the QND aspect in such a high-fidelity read- 
out is not expected to be an issue in the measurement of 
the final state of a quantum algorithm. However, reduc- 
tion in the QND character is problematic for tasks such as 
measurement-based state preparation, quantum feedback 
control and quantum error correction. 

In summary, we have shown that for a qubit with M > 
2 levels dispersively coupled to a resonator, the qubit- 
induced nonlinearity of the resonator depends in a nontriv- 
ial way on the qubit state. This can be exploited to increase 
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FIG. 5: Qubit relaxation (full black lines) and leakage rate 
(dashed blue lines) due to Purcell effect (a) and bare qubit de- 
cay (b) as a function of the measurement power. The dotted red 
line in (a) is the expected value for 7k at low power. The parame- 
ters are the same as in Fig[T] 



the SNR ratio in a QND measurement at low-photon num- 
ber and captures the essential aspects of the high-fidelity 
non-QND measurement recently reported 12]. 
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Note added. - Theoretical modeling of the high-fidelity 
readout has also been reported by Bishop et al. (following 
Letter) iH. 
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